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[i] The Geoelectrodynamics and Electro-Optical Detection of Electron and Suprathermal 
Ion Currents (GEODESIC) sounding rocket encountered more than 100 filamentary density 
cavities associated with enhanced plasma waves at ELF (<3 kHz) and VLF (3-10 kHz) 
frequencies and at altitudes of 800-990 km during an auroral substorm. These cavities were 
similar in size (~20 m diameter in most cases) to so-called lower-hybrid cavities (LHCs) 
observed by previous sounding rockets and satellites; however, in contrast, many of the 
GEODESIC cavities exhibited up to tenfold enhancements in magnetic wave power 
throughout the VLF band. GEODESIC also observed enhancements of ELF and VLF 
electric fields both parallel and perpendicular to the geomagnetic field B 0 within cavities, 
though the VLF E field increases were often not as large proportionally as seen in the 
magnetic fields. This behavior is opposite to that predicted by previously published theories 
of LHCs based on passive scattering of externally incident auroral hiss. We argue that 
the GEODESIC cavities are active wave generation sites capable of radiating VLF waves 
into the surrounding plasma and producing VLF saucers, with energy supplied by cold, 
upward flowing electron beams composing the auroral return current. This interpretation is 
supported by the observation that the most intense waves, both inside and outside cavities, 
occurred in regions where energetic electron precipitation was largely inhibited or 
absent altogether. We suggest that the wave-enhanced cavities encountered by GEODESIC 
were qualitatively different from those observed by earlier spacecraft because of the 
fortuitous timing of the GEODESIC launch, which placed the payload at apogee within a 
sub storm-related return current during its most intense phase, lasting only a few minutes. 
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1. Introduction 

[ 2 ] Lower-hybrid cavities (LHC), or lower-hybrid solitary 
stmctures (LHSS), consist of localized enhancements of 
wave power near and above the lower-hybrid cutoff fre- 
quency. The waves reside in density depletions of several 
to tens of percent and are sometimes observed to heat ions to 
several eV [LaBelle et al, 1986; Kintner et al. , 1 992; Arnoldy 
et al , 1992; Vago et al , 1992]. These first reports were based 
on sounding rockets in the topside auroral ionosphere (see 
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review by Lynch et al [1999]). Later, LHSS were observed 
by the Freja satellite up to 1700 km [Eriksson et al , 1994; 
Dovner et al , 1994; Knudsen et al , 1998] and on the Cluster 
satellites throughout the magnetosphere [ Tjulin et al , 2003]. 
The satellite measurements show that LHSS occur in the 
plasmapause and plasmasphere as well as in the auroral 
ionosphere. Their formation appears to depend on the pres- 
ence of ambient VLF hiss of sufficient amplitude [ Lynch 
et al , 1999]. Hoymork et al [2000] used chord length dis- 
tributions to show that LHSS are cylindrical as opposed to 
planar. 

[ 3 ] Interferometric analysis of LHSS wave electric fields 
has demonstrated, on at least three different sounding rockets 
and on Cluster, that waves above (below) the local lower- 
hybrid frequency^// propagate azimuthally in a right- (left-) 
handed sense around the cavities [. Pingon et al , 1997; 
Bonnell et al , 1998; Tjulin et al , 2003]. This behavior is 
consistent with the predictions of a passive scattering model 
in which linear mode conversion of externally incident 
whistler mode hiss drives short-wavelength lower-hybrid 
waves inside cavities [ Seyler , 1994; Schuck et al , 1998; 
Hall et al , 2004]. Similar examples of whistler-driven 
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Figure 1. Payload configuration, showing the placement of the electric field spheres (numbered), 
Langmuir probe (LP), search coil antennas (MFE), Suprathermal Electron Imager (SII), Thermal Electron 
Imager (TEI), and the Energetic Electron and Ion Detectors (EED and EID). The payload was aligned with 
the geomagnetic field to within 5° throughout most of the flight. 


lower-hybrid waves have been observed on density gradients 
in space [Bell and Ngo, 1990, and references therein; 
McAdams et al ., 1998] and in laboratory plasmas [ Bamber 
et al , 1995; Rosenberg and Gekelman , 2001]. Search coil 
magnetometer antennas have been used to confirm that the 
magnetic component of VLF emissions found in LHSS 
cavities and on density gradients is usually either absent or 
small, in general agreement with passive scattering theory 
[ Vago et al , 1992; Hoymork et al , 2000; Reiniusson et al , 
2006]. 

[ 4 ] Schuck et al [2003] provide a comprehensive review 
of LHSS. At the time of writing of that review many basic 
aspects of LHSS were explained as outlined above; how- 
ever, the cause of the density depletion was a notable 
exception. Earlier attempts to explain the depletions in terms 
of pondermotive expulsion along B 0 [Chang, 1993; Shapiro 
et al , 1993] were abandoned when it became apparent that 
observed wave amplitudes were insufficient to support the 
observed depletions [Robinson et al , 1996], some as deep as 
tens of percent [Kintner et al , 1992]. For that reason, LHSS 
depletions were assumed to be “preexisting” [Robinson et al , 
1996; Schuck et al , 2003], meaning they are due to an 
unspecified mechanism unrelated to the observed lower- 
hybrid waves and pointing to some additional process that 
must precede LHC formation. Knudsen et al [2004] showed 
that density depletions of tens of percent are a natural con- 
sequence of ion heating to several eV that is localized on the 
ion gyroradius scale as observed on several sounding rockets, 
including Geoelectrodynamics and Electro-Optical Detection 
of Electron and Suprathermal Ion Currents (GEODESIC). 

[ 5 ] In this paper we report observations by the GEODESIC 
sounding rocket of wave-filled cavities very similar to LHCs, 
but with a strong magnetic field enhancement that is not 
compatible with the passive scattering model. Moreover, 
the wave electric fields show none of the mode structure 
reported by Pingon et al [1997] and Tjulin et al [2003], 
indicating that the GEODESIC events are qualitatively 
different than LHC/LHSS. After describing the GEODESIC 
experiment and data in sections 2 and 3, we explore and 


discuss possible reasons for and the significance of these 
differences in sections 4 and 5. 

2. Experiment 

[6] The GEODESIC sounding rocket was launched at 
09:19:10.3 UT (00:19:10.3 LT, 22:10 MLT) on 26 February 
2000 from Poker Flat Research Range in Alaska. A Black 
Brant XII vehicle carried the payload to an apogee of 99 1 km. 
The payload spin axis was aligned by an attitude control 
system to within a few degrees of — B 0 ; the spin period was 
5.4 s. Figure 1 shows the configuration of the payload. 

[ 7 ] GEODESIC carried electron and ion top hat spectro- 
meters (0.01-30 keV), a suprathermal ion imager measuring 
the range 0-20 eV [see Burchill et al , 2004; Knudsen et al , 
2003], a fixed bias cylindrical Langmuir probe in ion satu- 
ration, and a swept Langmuir probe to measure electron 
density and temperature. 

[8] Electric fields were measured from DC to 3.5 MHz 
including broadband time series data and onboard fast Fourier 
transform power spectra and burst memory “snapshot” col- 
lection. For the data shown here, broadband wave data were 
gathered in a pass band of 0.015-10 kHz (Nyquist) with 
perpendicular double probes measuring 8.0 and 8.8 m tip to 
tip, oriented in the plane perpendicular to B 0 . Dual sensors 
on each forward boom provided shorter (1 m) baselines, as 
shown in Figure 1. An additional pair of booms measuring 
6 m tip to tip was deployed 2.2 m aftward, allowing a mea- 
surement of electric field components parallel to B 0 . 

[ 9 ] GEODESIC carried two orthogonal induction coil 
antennas (search coils) sensitive to magnetic perturbations in 
the plane perpendicular to the payload spin axis, and there- 
fore perpendicular to B 0 to within 5°. Each was 18 cm in 
diameter with 20,000 turns and contained no magnetic core 
material. Outputs were sampled at 20,000 s _1 . All search coil 
data presented in this study have been compensated accord- 
ing to a semiempirical response function that accounts for 
measured values of coil resistance, distributed capacitance 
and self-inductance. This response (output voltage per unit 
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Figure 2. Overview of GEODESIC measurements as a function of time after launch at 09:19:10.3 UT 
(00:19:10.3 LT), altitude, and corrected geomagnetic latitude (CGLat). (a) B±, power spectral density 
(PSD) measured by a search coil magnetometer sensitive to fluctuations in the plane perpendicular to B 0 . 

(b) E ± , electric field PSD measured over the short (1 m) “1-5” baseline in the plane perpendicular to B 0 . 

(c) E || , electric field PSD measured parallel to B 0 over the 2.2 m “5-9” baseline, (d) e- downgoing, electron 
energy flux parallel to B 0 . (e) Ion density, as measured from a fixed bias cylindrical Langmuir probe. 


input magnetic field) shows a peak near 500 Hz and a 
monotonic falloff proportional to l/f 2 at higher frequencies. 
In order to prevent contamination from electric fields, the 
coils (including input terminal pins) were placed inside 
grounded conductive enclosures. 

[ 10 ] Postflight analysis of the search coil data revealed an 
anomalously high degree of correlation between the two 
orthogonal channels at frequencies above ~1 kHz. The 
source of this correlation (electronic, inductive) remains 
unclear. With no physical model of the unwanted coupling, 
directional field information is compromised. Errors in 
absolute amplitude due to this coupling are also possible. 
Nevertheless, we are able to show in section 3 that the coils 


clearly respond to magnetic fluctuations of geophysical 
interest. We are also able to validate the magnetic calibration 
using expected properties of electromagnetic whistler waves 
(see section 3.2). 

3. Data 

3.1. Overview 

[li] Figure 2 summarizes the plasma and wave environ- 
ment through which GEODESIC flew. Starting at Figure 2e, 
the plasma density during the flight remained high (above 
10 4 cm -3 ) even at the 991 km apogee. Precipitating electrons 
shown in Figure 2d reach 20 keV within a broad inverted- V 
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Kaktovik magnetic field 



Figure 3. Ground-based magnetometer H (magnetic north) 
component at Kaktovic, Alaska. The vertical dashed line 
indicates the time at which GEODESIC was magnetically 
conjugate. After Burchill et al [2004]. 

structure in the early portion of the flight. The inverted-V 
energy rises abruptly just after 300 s, then breaks into softer, 
energy-dispersed bursts at ~330 s. Measurements at other 
pitch angles (not shown) demonstrate that these bursts are 
highly collimated along the magnetic field direction. A sec- 
ond, weak inverted-V is followed by another period with 
drastically reduced and intermittent downward precipitation 
lasting from approximately 470 to 570 s, then by a third 
inverted-V and finally by softer, structured precipitation 
characteristic of the poleward boundary of the auroral zone in 
the northern portion of the flight. 

[ 12 ] Figure 3 shows that the period from ~300-600 s after 
launch is marked by a strong —750 nT excursion in the H 
component magnetic field as measured on the ground at 
Kaktovic, overflown by GEODESIC at ~400 s after launch. 
Thus, the observed structure in the electron precipitation is 
strongly affected by the dynamics of a substorm expansion 
and is at least partly temporal. 

[ 13 ] Moving up to the electric field data in Figures 2c 
and 2b, spectrograms of both the perpendicular electric 
field E\\ (from baseline 1-5) and parallel electric field E\\ 
(baseline 5-9) are dominated by intense whistler mode 
auroral hiss with a pronounced lower cutoff between 5 and 
7 kHz, easily recognizable as the lower-hybrid frequency^//. 
Emissions below f LH are much weaker, with the exception 
of the period 480-560 s, where intense ELF wave emissions 
appear. We note that the most intense wave emissions in both 
the ELF and VLF are found in a region nearly (though not 
completely) devoid of electron precipitation. Absence of 
precipitation and strong ELF wave emission are commonly 
associated with downward current (upgoing electron) regions 
as identified at higher altitudes [Paschmann et al , 2003]. 
Furthermore, while failure of the attitude camera precludes 
an accurate DC magnetic field baseline measurement for 
GEODESIC, a partial solution (not shown) reveals a system 
of currents structured on scales of a few kilometers and with 
both positive and negative excursions of tens of //A m -2 . 
Because of these two observations, we identify the region 
traversed between 480 and 560 s as a likely site for intense 
auroral return currents. 

[ 14 ] Figure 2a demonstrates a significant magnetic com- 
ponent in the VLF hiss as measured by the search coil mag- 
netometer. While there are weak magnetic emissions near 
and above f LH throughout most of the flight, again the most 
intense wave emissions coincide with regions nearly devoid 
of electron precipitation, as in the case of the electric fields. 


The magnetic emissions differ qualitatively from the electric 
fluctuations in that there is no pronounced cutoff at f LH , 
although the magnetic spectrum intensity does decrease with 
decreasing frequency. Theoretically, no cutoff is expected, 
since the lower-hybrid cutoff applies to (essentially) elec- 
trostatic lower-hybrid waves propagating nearly perpen- 
dicular to B 0 , and not to R-O-mode electromagnetic hiss. A 
second characteristic of the magnetic data is the near absence 
of any signature from the ELF waves in the downward cur- 
rent region (480-560 s). 

[ 15 ] We note that noise measured on the launch pad just 
prior to launch was subtracted from Figure 2a. This included 
interference lines at frequencies of 4000 and 4375 Hz. We 
now turn to a study of small-scale wave-filled cavities found 
within the most intense VLF wave regions shown in Figure 2. 

3.2. Wave Cavities 

[ 16 ] The spectrograms and time series plots in Figure 4 
span 2.0 s during the flight at an altitude of 980 km. A 
strong magnetic interference signal at 4.375 kHz has been 
suppressed in Figure 4a (top). Five distinct wave cavities are 
visible, four of which exhibit strong magnetic signatures in 
both the spectral and time series displays in Figure 4a. The 
perpendicular electric fields in Figure 4b show clear, broad- 
band responses above the lower-hybrid (LH) cutoff within 
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Figure 4. Two-second sequence during the GEODESIC 
flight showing (a) B±, perpendicular magnetic, (b) E±, per- 
pendicular electric (baseline 1-5), (c) parallel electric 
(baseline 5-9), and (d) ion density signatures of ELF/VLF 
wave cavities. (Note that the density axis does not extend to 
zero.) 
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Figure 5. Red lines show power spectra averaged over 
15 cavities of field quantities defined in Figure 2. Blue lines 
are calculated using samples taken outside but adjacent to 
cavities. Error bars represent RMS variation over the 15 sam- 
ple spectra. E\\/E± shows the ratio of parallel and perpen- 
dicular field amplitudes; cB ± /E ± is a refractive-index-like 
quantity. Black lines show the theoretical values for whistler 
mode waves propagating at 0° (solid), 80° (dashed), and 89° 
(dot-dashed) with respect to B 0 . 


cavities, and even stronger signatures in the ELF. In this 
view the LH cutoff is not very clear even outside cavities; this 
is because of a strong ELF background in this region. 

[ 17 ] In contrast, the ELF response is strongly suppressed 
below f LH in is 1 1 (baseline 5-9). This is expected due to high 
field-aligned conductivity and is a strong indication that the 
measurement is, in fact, responding to parallel fields. 

[is] The times series of E ± in Figure 4b shows that the 
enhanced fluctuations are superimposed on a perturbation 
with a timescale of the order of the cavity width (~20 ms). A 
high-pass filter at 1 6 Hz on the channel shown here (VLF 1 5) 
has modified the appearance of the features. Inspection 
of a long-baseline, DC-coupled channel (not shown) shows 
them to be similar to the bipolar structures reported in 


LHCs by Vago et al [1992], who identified them as 
“payload-generated waves.” Using a pair of subpayloads 
connected by a long tether on the OEDIPUS-C rocket, 
Knuds en et al. [1999] detected bipolar LHC signatures and 
argued they could be explained by a DC potential of order 
—0.5 V inside cavities. Regardless of their cause, instru- 
mental or geophysical, the bipolar signature could affect the 
data only up to a few times the characteristic frequency of 
the cavities in the payload frame, that is, of the order of 100- 
200 Hz. In the data shown in Figure 4, electric field intensi- 
ties extend smoothly from tens of hertz all the way to 10 kHz. 
Apart from the possibility of payload-generated waves at the 
lowest end of this range, the ELF waves in the GEODESIC 
data appear to be geophysical. 

[ 19 ] Figure 4e shows plasma densities derived from a fixed 
bias Langmuir probe biased to —5 V. The depletion depths of 
the cavities shown here are of order 10%. A study of these 
depletions and a mechanism that can generate them is given 
by Knudsen et al. [2004]. Cavity widths based on ions, 
plasma density, and electric fields were reported by Burchill 
et al [2004], who identified 106 events in total. 

3.3. Cavity Wave Spectra 

[ 20 ] Figure 5 shows average spectra calculated from 15 
cavities having the strongest magnetic field intensities. Each 
panel compares average spectra measured inside these cavi- 
ties (in red) with spectra taken adjacent to them, 0.1 s away 
(in blue). Figure 5a shows that the average magnetic field 
power spectral density increases monotonically with fre- 
quency, with no discernible feature near f LH . (Points near 
4.375 kHz suffered from strong interference and were dis- 
carded and replaced with interpolated values.) The spectrum 
from cavity interiors has the same shape as that taken outside, 
but is 5-10 times more intense. 

[ 21 ] Figures 5b and 5c show averaged spectra of E ± and E\\ 
from the 1-5 and 5-9 baselines, respectively. The spectrum 
of the ambient E ± outside cavities (in blue) is almost flat and 
exhibits only a slight depression just below the cutoff fre- 
quency (5 kHz). Inside cavities, the E± power spectral den- 
sity (PSD) increases tenfold below f LH , and threefold to 
fivefold above, with no signature of the f LH cutoff. Presum- 
ably this cutoff is obscured by the intense ELF waves 
extending up to f LH [see also Burchill et al , 2004, Figure 7]. 

[ 22 ] Unlike the case for Ej_, the average E\\ spectrum out- 
side cavities is strongly attenuated below f LH . Inside cavities 
the emissions are again even more intense below f LH than 
above, and the average spectrum inside is similar to that 
of E±, though somewhat reduced. This is illustrated in 
Figure 5d for “E\\IE±_” showing a nearly constant spectral 
ratio of between 0.5 and 0.8 within cavities and a ratio of 
unity outside cavities and above f LH . 

[ 23 ] Figure 5e plots the refractive-index-like ratio cB ± /E ± . 
Normally, we would wish to compare components of E ± and 
B ± that are perpendicular to each other as well as to B 0 . But 
as discussed in section 2, orthogonal magnetic field channels 
appear to be coupled. This fact leads to uncertainties both 
in the direction and absolute amplitude of the magnetic 
fluctuations. Nevertheless, using single-channel calibration 
values determined prior to launch, Figure 5 shows that the 
ratio cB ± /E ± falls within the theoretical refractive index for 
whistler mode hiss [Stix, 1992] propagating parallel to B 0 
(solid black line.) This provides an independent check that 
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the magnetic calibration is reasonable to within a factor 
of ~2. 

[ 24 ] Inside cavities and above 1 kHz, the ratio cB ± /E ± is 
nearly equal to or slightly higher than its value outside cav- 
ities, on average. This result is independent of the absolute 
accuracy of the field measurements and depends instead only 
on their linearity. 

[ 25 ] As shown by the dashed lines in Figure 5e, whistler 
mode waves propagating nearly perpendicular to B 0 have a 
much smaller refractive index than measured either inside or 
outside the cavities. Since any waves localized within cavi- 
ties must have a highly oblique wave vector, the whistler 
mode is not a viable candidate to explain the magnetic field 
enhancements associated with cavities. 

4. Interpretation 

[ 26 ] GEODESIC’s detections of large magnetic enhance- 
ments are in marked contrast to all earlier observations, 
raising the question of their origin. In this section we explore 
an interpretation in which the field enhancements seen by 
GEODESIC were the result not of coupling to externally 
incident waves but instead were generated within the cavities 
themselves. We propose that unstable auroral return currents 
produce electrostatic waves having electric field components 
both perpendicular and parallel to B 0 . Fluctuating parallel 
electric fields then modulate the parallel current, leading to 
strong magnetic fluctuations. 

[ 27 ] A full description of the cavity fields shown in 
Figure 4 would require not only frequency but also wave 
number spectra that are simply not possible to obtain from a 
rapidly moving spacecraft that traverses the highly structured 
cavities in only ~20 ms and that can introduce significant 
Doppler shifts. This situation severely limits any theoretical 
description since the perpendicular ion response to wave 
fields depends strongly on wavelength and also on ion 
composition. On the other hand, the parallel response, being 
dominated by electron currents, is much simpler, especially 
in the case of the very low (3 plasma relevant to GEODESIC. 
(/? is the ratio of particle to magnetic pressure.) In this 
section we exploit this fact to arrive at simple and testable 
predictions relating perpendicular magnetic and parallel elec- 
tric field fluctuations. 

4.1. Source of Fluctuating B 

[ 28 ] Assuming cylindrical symmetry, the Maxwell- 
Ampere relation can be used to relate the average field- 
aligned current y'y over a circular cross section of radius 
R to the average perpendicular (azimuthal) magnetic field 
around its perimeter: 




Moil 


1 dE\ 


(i) 


Referring to Figure 5 a, the magnetic fluctuations in the 5- 
10 kHz range have an RMS amplitude of the order of 0.7 nT 
[(10 -4 nT 2 Hz -1 ) x (5000 Hz)] _1/2 . Ignoring the displace- 
ment current on the right-hand side of equation (1), a mag- 
netic field of this magnitude will result from a fluctuating 
field-aligned current having an RMS amplitude of ~ 140 //A 
m -2 over a 20 m diameter cylinder. Such a current could 
result, for example, from 20% of the background electron 


population (with a density of 10 4 cm -3 ) being accelerated 
to only 0.5 eV in the direction parallel to B 0 . 

[ 29 ] Next we take into account the first-order dynamics 
of field-aligned electron motion by separating parallel elec- 
tron velocity and electron density into background plus 
fluctuating components: v^y = v^jj >0 + n e = n e $ + 6n e . 
Field-aligned electron motion responds to a fluctuating par- 
allel electric field 8E\\ according to the parallel momentum 
equation 


+ (2) 

The field-aligned current is carried primarily by electrons: 
j || = q e n e v e ,\\. The first-order fluctuating current is 


8j\ | = q e {n e , o<5v e; || + 8n e v e ^) 


(3) 


[ 30 ] The fluctuating density and velocity are related by 
the continuity equation, dn e /dt + d(n e v e ^)/d || = 0. Assuming 
exp(/utf — ik\\z) dependence (z is the field-aligned coordinate), 


k\\n e $8v e ^\ 

V ~ k\\ v e, ||,0 


(4) 


Finally, combining equations (l)-(4), 

cB ± _ inR ( fp e 

i-v.h.o/v^ii ) 2 

where f pe is the electron plasma frequency, v ph> \\ = uj/k\\, c is 
the speed of light, and R is the radius of the current filament. 
The last term in (5) is due to parallel displacement current, 
and in this form can be seen to be less important than parallel 
conduction current by a factor of f pe /f 2 = 10 4 at 10 kHz for 
a density of 10 4 cm -3 . 

[ 31 ] Using observed values for n e 0 and R (10 m), the ratio 
v e ,\\,o/Vph,\\ is the only free parameter. Assuming an electron 
beam moving at four times its thermal speed (v e> y >0 = 4v^ e ), 
parallel wavelength Ay becomes the adjustable parameter. 
Figure 6 compares measured values of \cB ± /8E\\\ with 
equation (5) for various values of Ay, and shows that mea- 
sured values are consistent with Ay ~ 40-80 m over most of 
the spectrum, with a possible second branch of solution for 
longer wavelengths that could apply at the higher frequen- 
cies. Although this range of Ay is not verified independently, 
such wavelengths are plausible. 

[ 32 ] An additional test of the plausibility of this model 
comes from recognizing that if the magnetic fluctuations seen 
by GEODESIC are indeed due to fluctuating field-aligned 
currents confined within cylindrical density depletions, one 
would expect magnetic field intensity to be zero in the cavity 
centers. GEODESIC did not fly directly through cavity cen- 
ters in general, and so one would not expect to find zero (or 
even minimum) 8B in all crossings. However, some of the 
events do indeed exhibit a local minimum in magnetic wave 
power near their centers. Figure 7 shows an example as seen 
in a magnetic field intensity profile and wavelet spectrum, 
along with electric field intensity profiles using four different 
baselines. The differences in the electric field profiles are due 
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parallel wavelength 



Figure 6. (top) Refractive-index-like quantity (equation 
(5)) based on perpendicular magnetic and parallel electric 
fields averaged over 15 spectra measured inside (red) and 
outside (blue) cavities. The solid black line plots equation (5) 
for cavity radius R = 10 m, an electron beam velocity 
4 times its thermal speed, a beam temperature of 0.1 eV, and 
parallel wavelength Ai = 70 m. The dotted black line is 
an approximate fit of equation (5) to the measured values 
obtained by adjusting Ay. (bottom) The resulting parallel 
wavelength spectrum. The dashed black line in Figure 6 (top) 
illustrates a second branch of solution for a parallel phase 
velocity exceeding the electron beam velocity, in this case 
for Ay = 1000 m. 


to the fact that they sampled different cross sections of the 
cavity. 

[ 33 ] As noted, the simple model presented here does not 
include particle motion or currents perpendicular to B 0 , and 
therefore cannot describe a normal mode of the plasma. 
In particular, the magnetic fields predicted by this model 
should not be interpreted as the magnetic signature of 
whistler mode waves. Instead, the model provides a testable 
relation between two observables in the near field of an 
oscillating current filament driven by electrostatic waves 
having an electric field component parallel to B 0 . Calculating 
the coupling of these oscillations to whistler waves in the 
surrounding plasma would require a more complete plasma 
model that is beyond the scope of this study. 


Table 1. Input Parameters for WHAMP Wave Growth 
Calculations 


Species 

Mass 

(H) 

Charge 

(f«0 

Density 

(cm- 3 ) 

Temp 

(eV) 

V\\ drift 

Ions (0 + ) 

15.9281 

+1 

14 x 10 3 

0.2 

0.0 

Electrons 

0.0 

-1 

7-12.6 x 10 3 

0.2 

0.0 

Electron beam 

0.0 

-1 

1.4-7 x 10 3 

0.05-0.5 

2. 0-7.0 


observation. The lower boundaries of density cavities are 
possible regions of such acceleration. Because field-aligned 
current must be conserved, the drop in density as electrons 
enter cavities must be compensated by increased upward 
velocity. 

[ 35 ] To test the feasibility of this wave generation mecha- 
nism for GEODESIC, we estimate electrostatic wave growth 
rates for GEODESIC conditions using the WHAMP dis- 
persion solver [ Ronnmark , 1985], which assumes infinite 
plane waves in a homogeneous plasma. These assumptions 
clearly do not hold for the GEODESIC cavities. However, 
the density cavities are of the order of a full perpendicular 
wavelength in diameter, and in the absence of a more rep- 
resentative model we believe that WHAMP ’s estimates are 
valid in an order-of-magnitude sense, and are useful for 
establishing trends. 


a) 


N 


GEODESIC 





521.540 521.550 521.560 521.570 521.580 

Time after launch (s) 


4.2. Source of Fluctuating E\ \ 

[ 34 ] James [1976] argued that saucer-shaped VLF wave 
signatures are driven by electrostatic waves arising from a 
cold-electron drift instability in downward current regions. 
We extend this scenario one step further to say that parallel 
electric fields associated with the electrostatic waves modu- 
late the current-carrying beam and also drive any ambient 
electrons along B, generating a magnetic signature. Such a 
mixture of ambient (nondrifting) and beam electrons will 
occur when beam electrons are accelerated below the point of 


Figure 7. (a) Wavelet spectrogram and integrated power 
profile of magnetic fluctuations B ± across a single cavity, 
showing a minimum near the center, (b) Similar displays 
for electric fields obtained from the four 1 m long baselines 
perpendicular to B 0 . “VLF” indicates that these electric field 
channels are AC-coupled with a filter cutoff near the Nyqvist 
frequency of 10 kHz. (c) Ion density from the fixed bias 
Langmuir probe. The dark strip in Figure 7a (top) is due 
to suppression of a strong interference signal originating on 
the payload. 
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WHAMP Growth Rate 



Figure 8. Wave growth rates calculated from WHAMP for 
a variety of electron beam densities and temperatures, assum- 
ing a background density of 10 4 cm -3 and a perpendicular 
wavelength of 1 0 m. 

[36] Table 1 shows assumed parameters of the ambient 
plasma and the electron beam. Figure 8 shows strong wave 
growth for beam velocities of only a few times the electron 
thermal speed, a condition that is presumably easy to meet 
within a strongly driven ionospheric return current, as it 
implies electron energies of the order of only 1 eV or even 
less. Figure 9 illustrates that the drift of cold, current-carrying 
electrons can reasonably lead to parallel electric fields 
that are comparable to or even larger than perpendicular 
fields for wave numbers and other parameters typical of the 
GEODESIC environment. This is consistent with Figure 5d, 
which shows that E ± >E y near the lower-hybrid frequency as 
expected, and the two components become comparable at 
higher frequencies. 

[ 37 ] The argument above offers an explanation of the 
source of the strong waves observed by GEODESIC, but 
not the cause of their localization. Figure 10 shows the 


Growth Rate vs. Beam Density 



Figure 10. Dependence of electrostatic wave growth rate 
on relative beam density for parameters given in Table 1, 
supporting the idea that depletions in background density can 
enhance wave growth. 

dependence of the linear wave growth rate as a function 
of beam-to-ambient electron density as determined by 
WHAMP. For the same amount of DC field- aligned current, 
growth rate rises sharply as background density decreases. 
This suggests that fields both parallel and perpendicular to 
B 0 will grow preferentially in regions of locally depressed 
density. As shown by Knudsen et al [2004], locally 
enhanced perpendicular electric fields can enlarge ion gyro- 
orbits through heating so as to decrease time-averaged 
plasma density in the higher-field region. Since this depleted 
density can enhance instability according to Figure 10, a 
feedback loop can result, favoring density-depleted struc- 
tures localized on the scale of the heated ion gyroradii, of 
order 10 m. A self-consistent analysis of this mechanism 
is beyond the scope of this study; however, the simple 
model shows that localized ion heating to a few eV as seen 
on GEODESIC [ Burchill et al , 2004] can easily produce 



Figure 9. Electrostatic wave growth rates (color scale) in a homogenous plasma for parameters given 
in Table 1. (left) Growth rat q f of 10 kHz waves showing strong growth of waves having A^ = 10 m, 
Ay ~ 100 m, with corresponding values of | EJE X \ = \E\\/E ± \ > 1. (right) Wave growth rate as a function of 
real frequency f r , A_l, and Ay. 
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localized density depletions of the order of 10% [Knudsen 
et al ., 2004]. 

5. Discussion and Conclusions 

[ 38 ] The GEODESIC sounding rocket measured strong 
magnetic field fluctuations inside localized density cavities 
in the auroral ionosphere. Previous searches found that the 
magnetic component is either absent or weak [ Vago et al , 
1992; H 0 ymork et al , 2000; Reiniusson et al , 2006]. 

[ 39 ] In considering why GEODESIC found such a 
difference, the most likely explanation has to do with 
the geophysical environment. Whereas previous rocket- 
measured LHCs were found in or near auroral arcs, most of 
the GEODESIC events were found in a region of strongly 
reduced electron precipitation and highly structured and 
intense field-aligned currents, that is, a region in which 
intense auroral return currents are likely to be found. The 
magnetogram shown in Figure 3 is a signature of a substorm- 
related auroral electrojet and is not a direct measure of a 
return current, but it does indicate the presence of very 
intense ionospheric currents in the brief period during which 
GEODESIC passed over Kaktovic, on the north coast of 
Alaska. 

[ 40 ] In addition to the strong VLF magnetic fields, a sec- 
ond indication that the GEODESIC wave cavities are quali- 
tatively different from those reported previously is the lack 
of the telltale counterrotating eigenmode structure seen in 
frequency-wave number spectra of electric fields mea- 
sured by the AMICIST sounding rocket [Pingon et al , 1997; 
Bonnell et al , 1998]. This was determined by processing 
fields from several of the clearest GEODESIC events using 
the same procedure and software that were developed for the 
AMICIST rocket. This absence of resonant mode structure 
is consistent with a strong driving source within cavities. 

[ 41 ] Yet another unique aspect of the GEODESIC mea- 
surements is the high plasma density (~10 4 cm -3 ) near 
apogee, 1 to 2 orders of magnitude higher than is typically 
observed in the nighttime topside ionosphere. In the inter- 
pretation that electrostatic waves drive ambient electrons 
parallel to B 0 , the amplitude of the resulting magnetic fluc- 
tuations is proportional to electron density. It is conceivable 
that previous missions flew through wave generation regions 
having lower plasma density and therefore undetectable 
magnetic enhancements. 

[ 42 ] Strong parallel electric fields have been reported from 
the FAST satellite at higher altitudes [ Ergun et al , 2001], 
also in return current regions. High time resolution waveform 
snapshots on FAST show that both parallel and perpendicular 
electric field enhancements are due to discrete, charged “phase- 
space holes” traveling with the electron beam. GEODESIC 
waveform measurements were not able to resolve such 
structures; however, it is possible that the low-frequency 
E field signatures seen by GEODESIC are also the result 
of phase-space holes, which might be interpreted as the 
nonlinearly saturated state of the linearly unstable wave 
modes characterized by Figures 8 and 9. 

5.1. Measurement Anomalies and Uncertainties 

[ 43 ] As described in section 2, GEODESIC carried two 
orthogonal search coil antennas in order to determine the 
polarization of the magnetic waves in the plane perpendicular 


to B 0 . However, the two components exhibit an anomalously 
high degree of correlation at frequencies above 1 kHz, 
indicating an unwanted coupling of unknown origin and 
rendering polarization measurements unreliable. Coupling 
between the two channels will also affect amplitude mea- 
surements, though a simple model (namely, a cross-sensor 
coupling factor of less than unity) suggests errors of a fac- 
tor of two or less. This level of uncertainty is adequate to 
support the conclusions of this study. 

[ 44 ] The parallel electric field measurements summarized 
in Figures 2 and 4 are new and unusual. While parallel 
electric fields in return current regions are both predicted 
by theory and observed, at least at higher (FAST) altitudes, 
Figure 2 shows significant wave power in E\\ throughout 
the flight. Pickup of perpendicular fields due to imper- 
fect alignment with B is not a viable explanation since 
GEODESIC was aligned to within 5° throughout most of the 
flight, and the parallel component is comparable to and often 
stronger than E±, at least above the lower-hybrid cutoff fre- 
quency. Also, as illustrated in Figure 5, the shape of the 
perpendicular and parallel spectra are markedly different, 
meaning one channel cannot be explained simply through 
cross coupling with another. Improper boom deployment is 
also unlikely since perpendicular electric fields measured by 
the forward and aft boom sets agree well with each other. 
Finally, we note that our conclusions depend only on E\\ 
enhancements within cavities embedded in return current 
regions, and that an explanation of all aspects of the mea- 
surements shown in Figure 2 is beyond the scope of this 
study. 

[ 45 ] Regardless of the above issues and uncertainties, 
confidence in both B ± and E\\ measurements is enhanced 
considerably by the fact that together they are consistent with 
each other through a simple physical model of field-aligned 
current-carrying electrons driven by E\\ within cavities. 

5.2. Relation to VLF Saucers 

[ 46 ] VLF saucers are V-shaped features in frequency-time 
spectrograms observed by satellites over a wide range of 
altitudes [ Gurnett and Frank , 1972; Ldnnqvist et al, 1993]. 
Using ray tracing, James [1976, p. 514] inferred typical 
source dimensions of “0.5 km horizontally by no more than 
10 km vertically.” The total radiant power is of the order of 
1 0 mW, and the source spectrum is flat within a few decibels 
from 1 00 Hz to 20 kHz. James argued that saucers are driven 
by superthermal electrons with energies less than 5 eV which 
destabilize electrostatic waves as assumed in section 4 of 
this paper. 

[ 47 ] Using the FAST satellite, Ergun et al [2003] estab- 
lished that most but not all saucer sources are embedded 
in regions of downward current/upward electrons. In that 
study the upgoing (downgoing) electron variants had ener- 
gies above (below) 10 eV. They demonstrated that saucers 
are, in fact, composed of many individual “point sources” 
distributed in altitude and more confined in latitude than in 
longitude. Ergun et al [2001] also reported a strong corre- 
lation between the occurrence of saucers and electron phase 
space holes. 

[ 48 ] The altitude of the GEODESIC wave events, their 
location within a region in which auroral return currents are 
likely to occur, their highly localized nature, their relatively 
flat field spectra and the evidence that they are active sources 
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of VLF wave power are all consistent with these earlier 
reports of saucer source properties. Furthermore, the VLF 
power available within the GEODESIC cavities is sufficient 
to account for the James [1976] estimate of 10 mW for total 
radiant power. This can be seen by taking the product of a 
typical parallel electric field of 20 mV m -1 (0.1 (mV m -1 ) 2 
Hz -1 over a bandwidth of 5 kHz, see Figure 5) and typical 
field-aligned current fluctuations of 140 /a A m~ 2 (see 
section 4), resulting in over 800 mW of power available 
within a 20 m diameter, kilometer long cavity. A radiation 
efficiency of even 1% results in a total radiant power of 
the order of 10 mW as reported by James [1976]. 

[ 49 ] It would be preferable to observe directly a saucer 
resulting from the GEODESIC events; however, this would 
require a second, simultaneous observation far above the 
source. On the other hand, it is conceivable that some dis- 
persion might be visible from waves originating from within 
a single cavity a few km below GEODESIC’s trajectory; 
this was investigated in detail by Kabirzadeh [2010], and 
while highly localized saucer-shaped signatures were visible 
around events such as those shown in Figure 4, data pro- 
cessing artifacts could not be ruled out as an alternative 
explanation, and the results were inconclusive. 

[ 50 ] The results of this paper are as follows: 

1. GEODESIC sounding rocket observations show sev- 
eral-fold increases in VLF magnetic field amplitude localized 
within filamentary density depletions of the order of 1-10%. 

2. VLF electric field amplitudes both perpendicular and 
parallel to B 0 are also enhanced within the cavities, though on 
average the fractional increase is not as large as for magnetic 
fields. This observation excludes an interpretation in terms of 
lower-hybrid waves driven by linear coupling to externally 
incident whistler mode hiss as developed in previous studies 
ofLHSS. 

3. Observed magnetic field amplitudes perpendicular to B 0 
are consistent with an interpretation in which filamentary 
field-aligned currents were carried by electrons whose 
velocities were modulated by the observed parallel electric 
fields. 

4. VLF and ELF electric fluctuations are consistent with 
an interpretation in terms of filaments of unstable auroral 
return current carried by cold ionospheric electrons. 

5. The magnetic wave-filled cavities are found at altitudes 
between 800 and 990 km, within a region where downward 
electron precipitation is strongly reduced, and at the time of 
the most intense few minutes of an auroral substorm expan- 
sion as indicated by a —750 nT deflection of the magnetic H 
component as measured on the ground. 

6. A wavelet analysis identical to that used on data from 
previous sounding rockets and satellites (not shown here) 
shows no evidence for counterrotating eigenmodes, and 
therefore for passive scattering, supporting the conclusion 
that GEODESIC cavity waves were not passively scattered 
from waves originating outside cavities but instead were 
driven locally. 

7. The altitude, degree of structure, likely association with 
regions of auroral return current, spectral shape, and total 
power of the wave-filled cavities are consistent with proper- 
ties of VLF saucer source regions reported previously 
[James, 1976; Ergun et al. 2003]. 

[ 51 ] In summary, our results support a model of localized 
ELF/VLF wave cavities in which intense auroral return 


currents become unstable to electrostatic wave growth in 
the ELF/VLF frequency range when thermal electron flow 
becomes supersonic; these waves can energize ions and 
cause density depletions localized on scales comparable to 
heated ion gyro radii according to the mechanism reported 
by Knudsen et al [2004]. The depleted density destabilizes 
the waves further, expelling still more ions and leading to 
stable filamentary cavities. Fluctuating electric fields parallel 
to the geomagnetic field then drive electron currents parallel 
to B 0 ; these currents in turn generate magnetic oscillations as 
observed by GEODESIC. The fluctuating current filaments 
can radiate waves which disperse over large distances to 
form VLF saucers. When the intensity of the auroral return 
current diminishes, density cavities formed during the actively 
driven phase can persist and act as resonators for lower- 
hybrid waves excited by externally incident auroral hiss 
emissions, resulting in the passive lower-hybrid wave-filled 
cavities observed and described in previous studies. 
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